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Abstract. In this chapter we first calculate the post-Newtonian gravito-electric secular rate of the mean anomaly of 
a test particle freely orbiting a spherically symmetric central mass. Then, we propose a novel approach to suitably 
combine the presently available planetary ranging data to Mercury, Venus and Mars in order to determine, 
simultaneously and independently of each other, the Sun's quadrupole mass moment J20 and the post-Newtonian 
gravito-electric secular advances of the perihelion and the mean anomaly. This would also allow to obtain the 
PPN parameters 7 and j3 independently. We propose to analyze the time series of three linear combinations of the 
observational residuals of the rates of the nodes f2, the longitudes of perihelia -k and mean anomalies M of Mercury, 
Venus and Mars built up in order to account for the secular precessions induced by the solar oblateness and the 
post-Newtonian gravito-electric forces. According to the present-day EPM2000 and DE405 ephemerides accuracy, 
the obtainable precision would be of the order of 10 _4 -10 -5 for the PPN parameters and, more interestingly, 
of 10~ 9 for J20. Second, we investigate the possibility of measuring the post-Newtonian Lense-Thirring secular 
effect induced by the gravito-magnetic field of the Sun on the planetary orbital motion by analyzing two linear 
combinations of the residuals of the nodes f2 of Mercury, Venus and Mars. The proposed combinations are, by 
construction, unaffected by the quadrupole mass moment J2 of the Sun's gravitational potential. Moreover, they 
are insensitive also to the post-Newtonian gravito-electric field of the Sun because it affects only the perihelia u> 
and the mean anomalies M. The obtainable observational accuracy in the proposed measurement would be of 
the order of 19% and 36%. However, it must be pointed out that a major source of systematic error could be 
represented by the residual mismodelled classical JV— body secular precessions. The second proposed combination 
is insensitive just to such aliasing effect. 

Key words. Relativity - Gravitation- Celestial Mechanics- Sun: fundamental parameters - Planets and satellites: 
general- Methods: miscellaneous 



1. Introduction 

In this chapter we wish to investigate the possibility of performing some tests of Newtonian and post-Newtonian 
gravity in the Solar System arena extending to it some approaches which, in recent years, have proven to be fruitful 
in certain Earth- artificial satellites systems. More precisely, wc will deal with certain classical and post-Newtonian 
gravito-electromagnetic secular effects affecting the longitudes of the ascending nodes tt, the longitudes of the perihelia 
7T and the mean anomalies M. of the Solar System's planets. The proposed tests could be done in a relatively simple and 
cheap way (with respect to new missions involving the construction and launch of one or more spacecrafts); moreover, 
they would also benefit of the improvements in the planetary ephemerides coming from future planned space missions 
like Messenger and BepiColombo. The obtainable accuracies are comparable with those of future proposed or planned 
expensive and complex interplanetary missions. In the case of the gravito-magnetic field, it would be the first attempt 
to measure its effect on the planetary orbital motion in the gravitational field of the Sun. 

Historically, one of the first tests of the Einstein General Theory of Relativity (GTR) was the successful explanation 
of the anomalous secular gravito-electric perihelion advance of Mercury in the gravitational field of Sun (Einstein 1915) . 
On the contrary, another prediction of the Einstein's theory of gravitation for the motion of test particles, i.e. the 
much smaller gravito-magnetic Lense-Thirring effect (Lense & Thirring 1918), is still awaiting for an undisputable 
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Table 1. Orbital parameters of Mercury, Venus and Mars (http://nssdc.gsfc.nasa.gov/planetary/factsheet/). 
For the Astronomical Unit (A.U.) we use the value 1 A.U.=149597870691 m of the DE405 (Standish 1998) and 
EPM2000 ephemerides (Pitjeva 2001a; 2001b). The angle e refers to the inclination to the ecliptic. 



Planet a (A.U.) e (°) e 

Mercury 0.38709893 7.00487 0.20563069 

Venus 0.72333199 3.39471 0.00677323 

Mars 1.52366231 1.85061 0.09341233 



and direct observational check. As we will see, these post-Newtonian features of the planetary motion are strictly 
connected with the problem of the quadrupole mass moment J20 of Sun. In (Pireaux and Rozelot 2003) a theoretical 
range J20 = (2± 0.4) x 10~ 7 is admitted for the Sun's oblateness. For this topic and the interplay between this effect 
and the gravito-electric Einstein precession see the recent review (Pireaux & Rozelot 2003) and the references therein 
and (Ciufolini & Wheeler 1995). Basically, the point is the following. In regard to the secular orbital motions, to 
which we are interested here, the Sun's mass quadrupole moment induces classical effects which have, qualitatively, 
the same temporal signature of the relativistic ones. This means that they could corrupt the recovery of the genuine 
post-Newtonian features of motion because they could not be removed from the time series without removing the post- 
Newtonian signal of interest as well. All depends on the precision with which J20 is known: should the mismodclling 
in it induce classical residual precessions larger than the relativistic ones, there would be no hope to get a reliable 
test of relativistic gravity. As we will see later, the same problems could come from the N— body secular precessions. 
In this chapter we propose to disentangle such effects by measuring them in an independent way. More precisely, we 
propose to extend a certain approach used in Earth artificial satellite motion analysis to the interplanetary arena 
in order to single out just certain post-Newtonian (and Newtonian) orbital motion features independently of the 
Parameterized-Post-Newtonian (PPN) framework (Will 1993) which is usually employed in testing competing metric 
theories of gravity. Instead, the current approach consists of testing the post-Newtonian equations of motion as a whole 
in terms of the PPN parameters which are determined from multi-parameters fits together with other astrodynamical 
quantities. 

In Tabled the orbital parameters of the inner planets of the Solar System are reported. 



2. The post-Newtonian effects 

2.1. The gravito-electric effects 

In the framework of the standard PPN formalism the post-Newtonian gravito-electric acceleration induced by the 
Schwarzschild-like part of the spacetime metric and experienced by a test body freely falling around a static, spherically 
symmetric central mass M is (Soffcl 1989) 



GM 

a GE = -5—7 



2(7 + /3) j(v-v) 



r + 2(l+7)(r-t))i) , (1) 



where G is the Newtonian gravitational constant, c is the speed of light in vacuum, r and v are the position and velocity 
vectors, respectively, of the test body, 7 and (3 are the standard Eddington-Robertson-Schiff PPN parameters (7 = 
j3 = 1 in GTR; in this case eq.JQ) reduces to the expression of (Mashhoon et al. 2001)). Note that r here is the standard 
isotropic radial coordinate, not to be confused with the Schwarzschild radial coordinate r — r[l + GM / (2c 2 r)] 2 . In 
the usual orbital data reductions it is just r which is employed. 

By considering it as a small perturbation to the Newtonian monopole acceleration it is possible to work out its 
effect on the orbital motion of a test body with the standard perturbative techniques. The secular rate of the argument 
of pericentre is given by the well known formula 

ZnGM 

^ GE - Afl-e 2 )^' ( } 

where 

2 + 27-/3 ... 

^GE = ^ , (3) 



a and e are the semimajor axis and the eccentricity, respectively, of the test particle's orbit and n = ^/GM/a 3 is the 
(unperturbed) Keplerian mean motion. The orbital period of an unperturbed, two-body Keplerian ellipse is P = 2it /n. 
Note that eq.lfSJl is an exact result valid to all order in e. 



L. Iorio: The solar gravito-magnetic field and its possible measurement 



3 



As we will show, also the mean anomaly M. is affected by a post-Newtonian gravito-electric secular rate. The Gauss 
perturbative equation for M. is given by 

dM 2 A fr\ r / dui dfl\ , A . 



dt na \aJ \dt dt 

where f2 and i are the longitude of the ascending node and the inclination of the test particle's orbit to the equator 
of the central mass, respectively, and Ar is the radial component of the perturbing acceleration. In order to obtain 
the post-Newtonian gravito-electric secular rate of the mean anomaly one must consider the radial component of 
cq.£E}, insert it in the second term of the right-hand-side of eq.J3J, evaluate it on the unperturbed Keplerian ellipse, 
characterized by r = a(l — e 2 )/(l + e cos /) and by the radial and along-track components of the velocity vector which 
are vr — nae sin //Vl — e 2 , vt = na(l + ecos/)/Vl — e 2 , respectively (/ is the orbiter's true anomaly), multiplying 
it by 

dt = (1-e^df 

P 27r(l + ecos/) 2 ' [ ' 

and integrating over one orbital revolution, i.e. from to 2~k. It turns out that the post-Newtonian gravito-electric 
secular rate of the mean anomaly is, to order 0(e 2 ) 



nGM 
M GB 5 — 



We can define 



(2 + 47 + 3/5) ( 1 + - ) + (2 + 7 )e< 



(6) 



2 + 4 7 + 3/3 

MGE = 9 ' (7) 

It turns out that the second term of eq.jHJ) induces for Mercury an additional shift of 1.789 arcseconds/century (" 
cy _1 in the following), while for the other planets it is of the order of 10 -2 -10~ 4 " cy _1 . Thus, for them the secular 
rate of the mean anomaly can be written as 

9nGM 

M GE — ? ==MGE- (8) 

rayl — e 

However, as we will see later, the mean anomaly of Mercury will not be used in the combined residuals strategy outlined 
in the following. For Mercury, the post-Newtonian gravito-electric effect induced by the Sun on ui and A4, according to 
GTR, amounts to 42.980 " cy _1 and -127.949 " cy _1 , respectively. The most accurate estimate of the gravito-electric 
perihelion advance seems to be that obtained for Mercury by E.M. Standish in 2000 with the DE405 ephemerides 
(Standish 1998) and reported in (Pireaux and Rozelot 2003). He averaged the Mercury's perihelion evolution over 
two centuries by using the DE405 ephemerides with and without the post-Newtonian accelerations. Standish included 
in the force models also the solar oblateness with J2© = 2 x 10 -7 , so that the residuals for Mercury accounted for 
the post-Newtonian effects only; the determined shift was 42.98 ± 0.0023 " cy -1 . The same approach for the mean 
anomaly (Standish E M private communication) has yielded to —130.003 + 0.0027 " cy^ 1 (See also Table0]later). Note 
that the quoted uncertainty do not come from direct observational errors. They depend on the fact that in the force 
models used in the numerical propagation many astrodynamical parameters occur (masses of planets, asteroids, etc.); 
their numerical values come from multiparameter fits of real data and, consequently, are affected by observational 
errors. Such numerical tests say nothing about if GTR is correct or not; they just give an idea of what would be the 
obtainable accuracy set up by our knowledge of the Solar System arena if the Einstein theory of gravitation would be 
true. 

In regard to the possibility of constructing time series of planetary mean anomalies, it must be noted that a certain 
effort would be required. Indeed, so far this orbital elements has never been utilized, so that the partials, e.g., should 
be computed (Pitjeva, private communication 2004). 



2.2. The gravito-magnetic Lense-Thirring effect 

Another post-Newtonian secular precession which affects not only u> but also f2 is the Lense-Thirring effect (Lcnsc & 
Thirring 1918) induced by the proper angular momentum J of the central body 

— 6GJcosi 

^ LT = c 2 3 ( l_ e 2 ) 3/ a ^T, (9) 
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where 

MLT = (11) 

By assuming for the Sun J Q = 1.9 x 10 48 g cm 2 s _1 (Pijpcrs 2003), the gravito-magnetic effect on, e.g., Mercury's 
perihelion is of the order of 10~ 3 " cy _1 (De Sitter 1916). Such small value is at the edge of the present sensitivity 
(0.002 " cy _1 ) in determining the Mercury perihelion shift from the ephcmeridcs (cfr. the results obtained by Standish). 
Moreover, it should be considered that, even if future improvements of the obtainable observational sensitivity in inter- 
planetary ranging allowed to consider the possibility of measuring the Lcnsc-Thirring perihelion advance of Mercury, 
the impact of the systematic errors due to the uncertainties in the solar oblateness would severely limit the realistic 
accuracy obtainable in such demanding measurement. Indeed, by assuming an uncertainty of oj 20 = 0.4 x 10 -7 , the 
error in the secular precession induced by the Sun's quadrupole moment would amount to 245% of the Lense-Thrring 
shift of 7r Merc (see below for the definition of n). Another important source of systematic error in the measurement of 
such a tiny effect would be represented by the N— body classical secular precessions which are of the order of 10 2 — 10 3 
" cy _1 (see on the WEB http://ssd.jpl.nasa.gov/elem_planets.html). Indeed, their residual mismodelled part 
could severely bias the recovery of the Lense-Thirring effect. In this chapter we will propose a possible approach to 
overcome these problems. It must be noted that, if, on the other hand, we look at a Lense-Thirring test as a way to 
measure the Sun's angular momentum by assuming the validity of GTR, a measurement of the solar gravito-magnetic 
field would have a significance only if the obtainable accuracy was better than 10%; indeed, among other things, 
the present-day uncertainty in the Sun's angular momentum J Q , which could be measured from the Lense-Thirring 
precessions, is just of the order of 10% (Xu & Ni 1997) in various solar models (Paterno et al. 1996; Elsworth et al. 
1995) or even less in the framework of asteroseismogyrometry (Pijpers 2003). 

At present, the only performed attempts to explicitly extract the Lense-Thirring signature from the data of orbiting 
masses in the Solar System are due to Ciufolini and coworkers who analyzed the laser-ranged data of the orbits of 
the existing LAGEOS and LAGEOS II Earth artificial satellites (Ciufolini et al. 1998). A 20-30% precision level in 
measuring the terrestrial gravito-magnetic field is claimed, but other scientists judge these evaluations too optimistic 
and propose different error budget (Rics et al. 2003). In April 2004 the GP-B spacecraft has been launched. It will carry 
out a very complex and challenging mission which should be able to measure another gravito-magnetic effect, i.e. the 
precession of the spins of four superconducting gyroscopes (Schiff 1960) carried on board at a claimed accuracy of 1% 
or better (Everitt et al. 2001). Almost twenty years ago it was proposed to launch a third LAGEOS-like satellite— the 
LAGEOS III/LARES— and to analyze the time series of the sum of the residuals of the nodes of LAGEOS and LARES 
(Ciufolini 1986) or some other combinations of residuals of the nodes and the perigees of LARES and both the existing 
LAGEOS satellites (Iorio et al. 2002a). The obtainable accuracy would probably be of the order of 1% too. Mainly 
funding problems have prevented, up to now, from implementing such relatively easy and cheap mission. Recently, 
the possibility of measuring the Lense-Thirring precessions by means of the Relativity-dedicated OPTIS spacecraft, 
which could be launched in the same orbital configuration of LARES, has been considered (Iorio et al., 2004). The 
recently proposed LATOR (Turyshev et al. 2004) and ASTROD (Ni et al. 2004) missions would be sensitive to the 
gravito-magnetic part of the bending of light rays and time delay in the gravitational field of the Sun. Finally, it must 
be noted that, according to K. Nordvedt Jr. (Nordvcdt 2003), the multidecadal analysis of the Moon's motion with 
the Lunar Laser Ranging (LLR) technique strongly supports the existence of the gravito-magnetic force 1 as predicted 
by GTR, although in an indirect way. "It often has been claimed that the presence of gravitomagnetism within the 
total gravitational interaction has not been experimentally confirmed and measured. Indeed, different experiments 
have been under development to explicitly observe the effects of this historically interesting prediction of general 
relativity. But this gravito-magnetic acceleration already plays a large role in producing the final shape of the lunar 
orbit, albeit in conjunction with the rest of the total equation of motion; the precision fit of the LLR data indicates 
that gravitomagnetism's presence and specific strength in the equation of motion can hardly be in doubt. [...] It would 
be impossible to understand this fit of the LLR data without the participation of the gravito-magnetic interaction in 
the underlying model, and with strength very close to that provided by general relativity, 7 = 1". 

3. The solar oblateness 

The solar quadrupole mass moment J20 is an important astrophysical parameter whose precise knowledge could 
yield many information about the inner structure and dynamics of our star. A reliable evaluation of J 2 © still faces 
some controversy: on one side, the theoretical values strongly depend on the solar model used, whereas accurate 

1 According to Nordvedt (2003), the Earth-Moon range is affected by long-periodic harmonic perturbations of gravito- 
magnetic origin whose amplitudes are of the order of 5 m and the periods are monthly and semi-monthly. The amplitudes of 
the lunar motion at both these periods are determined to better than half a centimeter precision in the total orbital fit to the 
LLR data. 
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Table 2. Post-Newtonian gravito-electric precessions and coefficients of Newtonian precessions of the node, the peri- 
helion and the mean anomaly for Mercury, Venus and Mars in " cy _1 . The values 7i"GE) and -Mqe are calculated with 
GTR. The coefficients f2 2 , tt. 2 and M 2 are <9(f2j 2Q )/<9(J 2Q ), 9(7rj 20 )/9(J 2Q ) and d(Mj 2Q )/d(J 2 o), respectively. For 
fl and u> they can be found in (Iorio 2003). In order to have the precessions they must be multiplied by J 2 ©. Note 
that the result for the mean anomaly of Mercury accounts for the correction of order C(e 2 ) which, instead, can be 
neglected for the other planets. 



Planet J7ge itge Mge ^.2 ^.2 M.2 

Mercury 42.981 -127.949 -126878.626 126404.437 123703.132 

Venus 8.624 -25.874 -13068.273 13056.803 13056.504 

Mars 1.351 -4.035 -980.609 980.353 976.067 



measurements are very difficult to obtain from observations. For all this matter see the recent review (Pireaux & 
Rozelot 2003) and (Rozelot et al. 2004) . From an observational point of view, J 2 © is not directly accessible. In this 
context a dynamical determination of </ 2 ©, analyzing, e.g., the orbits of the inner planets of the Solar System, is 
interesting because it might be compared with those derived from solar model dependent values of the oblateness. 
However, it is not simple to reach this goal because of the interplay between the effects of the solar quadrupole moment 
with those induced by the post-Newtonian gravito-electromagnetic forces. Instead of the trajectory of planets, it would 
be possible to infer J 2 q from accurate tracking of some drag-free spacecraft orbiting within a few radii of the solar 
center. This will be the approach followed by, e.g. the BepiColombo mission (see Section l4~2l . Alternatively, the Sun's 
quadrupole mass moment can be inferred from in-orbit measurement of solar properties, like the SOHO-MDI space- 
based observations (Armstrong & Kuhn 1999), or from Earth-based observations like those realized, e.g., with the 
scanning heliometer of the Pic du Midi Observatory (Rozelot et al. 2004). 



3.1. The classical precessions induced by the solar oblateness 

For a given planet of the Solar System orbiting the Sun, apart from the classical effects induced by the precession 
of the equinoxes and by the other planets and major asteroids which are routinely accounted for in the ephemerides 
computations (Pitjeva 2001a; 2001b), the oblateness of Sun induces also secular precessions on f2, n = uj + Q cosi and 
M given by 

• 3 nJ 20 f Rq\ 2 , . 



J20 2 (1 - e 2 ) 2 



a 



3 nJ 2Q (Rq\ {3 2 . 



4 (1 -e 2 ) 3 / 2 



a 



The explicit expressions for the secular precessions of SI and u> up to degree I = 20 can be found in (Iorio 2003). For 
Mercury 2 the Newtonian precessions due to Sun oblateness, with J 2 q = 2 x 10~ 7 , are of the order of 10~ 2 " cy _1 for 
fi, 7r and M. 



4. The interplay between the solar oblateness and the post-Newtonian precessions 

In Table|2]the relevant parameters of the classical and post-Newtonian gravito-electric secular precessions of the nodes, 
the perihelia and the mean anomalies of Mercury, Venus and Mars are reported. In all the relativistic tests performed 
up to now by analyzing the perihelia advances only of the inner planets of the Solar System with the radar ranging 
technique (Shapiro et al. 1972; 1976; Shapiro 1990) it has been impossible to disentangle the genuine post-Newtonian 
gravito-electric contribution of eq.© from the Newtonian precession of ea. l|13|) . Indeed, the observational residuals 
of 7T for a single planet, built up by suitably switching off the post-Newtonian 0(c~ 2 ) terms and the oblateness of 

2 As pointed out in (Milani et al., 2002), the angle i refers to the inclination between the planet's orbital plane and the fixed 
reference plane of the celestial reference frame; it is not the angle e between the planet's orbital plane and the ecliptic. It turns 
out that i ~ e/2. For Mercury e = 7.00487°. 
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Table 3. Post-Newtonian gravito-magnetic precessions and coefficients of Newtonian precessions of the node for 
Mercury, Venus and Mars in " cy -1 . The values Olt are calculated with GTR. The coefficients Ctj, are d(Clj e )/9(J^ Q ). 
They can be found, up to degree I = 20, in (Iorio 2003). In order to have the precessions they must be multiplied by 

JlQ- 



Planet 




£1.2 


n. 4 


^N — body 


Mercury 


0.001008 


-126878.626476 


52.774935 


-446.30 


Venus 


0.000144 


-13068.273031 


1.349709 


-996.89 


Mars 


0.000015 


-980.609460 


0.023554 


-1020.19 



Sun in the force models of the equations of motion in the orbital processors softwares, account entirely for the post- 
Newtonian and the solar oblateness 3 effects. This is a unsatisfactory situation, both if we are interested in testing 
post-Newtonian gravity and if we want to obtain a dynamical, model-independent measurement of J20. Indeed, it is, 
of course, impossible to constraint both the effects if only one perihelion rate is examined one at a time: in recovering 
one of the two effects we are forced to consider the other one as if it was known. Since the post-Newtonian gravito- 
electric effect is three orders of magnitude larger than that induced by solar oblateness, a determination of the latter 
by assuming the validity of GTR would be affected by a non negligible systematic error induced by the precision 
to which the post-Newtonian pericentre advance is known from other (more or less indirect and more or less biased 
by other aliasing effects) tests (Lunar Laser Ranging, binary pulsars periastron advance 4 ). The inverse situation is 
more favorable: indeed, if we are interested in the post-Newtonian gravito-electric effect the relative systematic error 
induced on its measurement by the precession due to the solar oblateness amounts to 5 x 10~ 4 even by assuming for 
the latter effect a 100% uncertainty 5 . 

4.1. The present-day approach to test post-Newtonian gravity 

At this point it may be interesting to clarify what is the current approach in testing post-Newtonian gravity from 
planetary data analysis followed by, e.g., the Jet Propulsion Laboratory (JPL). In the interplay between the real data 
and the equations of motions, which include also the post-Newtonian accelerations expressed in terms of the various 
PPN parameters, a set of astrodynamical parameters, among which there are also 7 and /?, are simultaneously and 
straightforwardly fitted and adjusted and a correlation matrix is also released. This means that the post-Newtonian 
equations of motion are globally tested as a whole in terms of, among other parameters, 7 and [3; no attention is paid 
to this or that particular feature of the post-Newtonian accelerations. The point is that the standard PPN formalism 
refers to the alternative theories of gravitation which are metric, i.e. based on a symmetric spacetime metric. But 
it is not proven that an alternative theory of gravitation must necessarily be a metric one. Moreover, even in the 
framework of the metric alternative theories, the PPN formalism based on 10 parameters is not sufficient to describe 
every conceivable metric theory of gravitation at the post-Newtonian order; it only describes those theories with a 
particularly simple post-Newtonian limit. One would, in principle, need an infinite set of new parameters to add to 
the standard ten parameter PPN formalism in order to describe the post-Newtonian approximation of any a priori 
conceivable metric theory of gravity (Ciufolini 1991; Ciufolini & Wheeler 1995). 

4.2. The possibilities opened by the future missions 

Concerning the possibility of disentangle the effects of the solar oblateness from those of the post-Newtonian gravito- 
electric force, it is stated that the future space mission BepiColombo 6 of the European Space Agency (ESA) will 

3 Also if the solar oblateness is included in the force models, the related uncertainty induces a corresponding systematic error 
in the recovered post-Newtonian effect. Fortunately, it is small; by assuming o"j 2Q = 0.4 x 10~ 7 , it amounts to 0.01%. 

4 Note that the binary pulsars periastron measurement should not be considered as a test of relativistic gravity because the 
masses of the binary system are not known (Stairs et al. 1998) . The Lunar Laser Ranging measurements do not allow to single 
out uniquely the gravitoelctric pericentre advance from the other post-Newtonian features of motion of the Earth-Moon system 
(Nordvcdt 2001). Recently, it has been proposed to measure the relativistic gravito-electric perigee advance of the terrestrial 
LAGEOS II satellite (Iorio et al. 2002b; Lucchesi 2003), but, up to now, the experiment has not yet been performed. 

5 Note that the impossibility of disentangle the gravito-electric and Lense-Thirring effects would not seriously affect the 
recovery of the gravitoelctric precession: indeed, the bias induced by the gravito-magnetic effect on the gravito-electric shift 
amounts to 0.004% only for Mercury. 

6 See on the WEB http://astro.estec.esa.nl/BepiColombo/. Present ESA plans are for a launch in 2010-2012. 
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Table 4. Present-day accuracy in determining the node, perihelion and mean anomaly secular rates of Mercury, Venus 
and Mars according to DE405 (Standish 1998) and EPM2000 (Pitjeva 2001a; 2001b) ephemerides. The figures, in " 
cy^ 1 , represent the formal, statistical errors. Realistic errors should be 10 times larger, at least. While the results by 
Standish come from the mathematical propagation of the nodes, the perihelia and the mean anomalies (Standish, E. 
M. private communication) evolution with and without post-Newtonian terms (with 7 = (5 = 1) and their average 
over a time span of two centuries, the results by Pitjeva are based on real data. The figure for Mercury has been 
obtained in 2001, while the other ones have been determined subsequently (Pitjeva, personal communication 2004). 



Planet 


a^JDE405) 


**caic(££405) 




cr* obs (£PM2000) 


Mercury 


0.000182 


0.0023 


0.0027 


0.0086 


Venus 


0.000006 


0.0414 


0.0414 


0.1037 


Mars 


0.000001 


0.0014 


0.0014 


0.0001 



provide us, among other things, with a dynamical, model-independent and relativity-independent measurement of J2© 
by measuring with high precision the nodal motion of Mercury (Milani et al. 2002; Pireaux & Rozelot 2003) which 
is not affected by the post-Newtonian gravito-electric force. The claimed accuracy would amount to <tj 20 = 2 x 10~ 9 
(Milani ct al. 2002). However, such evaluation refers to the formal, statistical obtainable uncertainty only. Indeed, 
the residuals of the Mercury's node would account, to a certain level of accuracy, for the Lense-Thirring precession 
as well. By considering such effect as totally unmodelled in the force models, its impact on the measurement of J2© 
would induce a 8 x 10~ 9 systematic error. The formal, statistical accuracy for 7 and f3 is evaluated to be of the order 
of 2 x 10~ 6 (Milani et al. 2002). Also the ESA astrometric mission GAIA 7 should measure, among other things, the 
solar quadrupole mass moment by analyzing the longitudes of the ascending nodes of many minor bodies of the Solar 
System. The obtainable accuracy for 7 is of the order of 10~ 5 -10~ 7 (Vecchiato et al. 2003). The ASTROD mission 
should be able to measure J20 with a claimed accuracy of the order of 10~ 8 or, perhaps, 10~ 9 — 10~ 10 (Ni et al. 
2004). The claimed obtainable accuracy for the PPN parameters is 4.6 x 10~ 7 for 7 and 4 x 10~ 7 for (3. Further 
improvements may push these limits down to 10~ 8 — 10~ 9 . The recently proposed LATOR mission should be able to 
measure, among other things, 7 to a 10~ 8 accuracy level and J20 to a 10 -8 level (Turyshev et al. 2004). However, 
LATOR, whose goal is an extremely accurate measurement of the light deflection at the Sun's limb, should not allow 
to disentangle the various post-Newtonian contributions (first and second order gravito-electric contributions, first 
order gravito-magnetic term, solar oblateness effect) to it. The solar orbit relativity test SORT (Melliti et al. 2002), 
which would combine a time-delay experiment with a light deflection test, should allow to reach a 10 -6 accuracy in 
measuring the PPN parameters 



5. The observational accuracy in planetary radar ranging 

Concerning the present-day accuracy of the planetary radar ranging, the radar itself is accurate well below the 100 
m level (Standish 2002). The problem, however, comes from the fact that the surfaces of the planets have large 
topographical variations. They are modeled in different ways. For Mercury, spherical harmonics and some closure 
analysis (comparing values when two different measurements reflect off from the same spot on the surface) have been 
done (Anderson ct al. 1996) in DE405. For Venus, a topographical model, which comes from (Pcttcngill et al. 1980), 
has been used. For Mars, closure points can be used. Closure points are pairs of days during which the observed points 
on the surface of Mars are nearly identical with respect to their longitudes and latitudes on Mars. Since the same 
topographical features are observed during each of the two days, the uncertainty introduced by the topography may 
be eliminated by subtracting the residuals of one day from the corresponding ones of the other day. The remaining 
difference is then due to only the cphemeris drift between the two days. The closure points for Mars have a priori 
uncertainties of about 100 m or less when the points are within 0.2 degrees of each other on the martian surface. Of 
course, for Mars, there is also the spacecraft ranging - far more accurate than the radar: Viking Landers (1976-82), 10 
m; Mars Global Surveyor and Odyssey (1999-2003), 2-3 m. However, correction for Mars topography is possible not 
only by using closure points (in this method no all observations may be used), but with help of modern hypsometric 
maps and by the representations of the global topography with an expansion of spherical functions (Pitjeva 2001b). 



7 See on the WEB http://astro.estec.esa.nl/GAIA/. Present ESA plans are for a launch in mid-2010. 
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6. The proposed approach for disentangle the solar oblateness and the post-Newtonian effects 

In regard to the possibility of measuring explicitly some post- Newtonian effects from the analysis of the orbital motion 
of proof masses in the gravitational field of a real (rotating) astronomical mass like the Earth, it must be pointed 
out that the main problems come from the aliasing effects induced by a host of classical orbital perturbations, of 
gravitational and non-gravitational origin, which unavoidably affect the motion of the probes along with GTR. In 
particular, the even zonal harmonics 3g of the multipolar expansion of the gravitational potential of the central mass 
induce secular classical precessions which, in many cases, tend to alias the gravito-electromagnetic ones of interest. 
Moreover, also the non-gravitational perturbations, to which the perigees of the LAGEOS-like satellites are particularly 
sensitive, are another important source of bias. The approach proposed by Ciufolini (1996) and Iorio (Iorio 2002; Iorio 
& Morea 2004) in the performed or proposed tests with LAGEOS and LAGEOS II consists of suitably designing linear 
combinations ^2 u «^^obs + S w j<^obs °^ or bital residuals which are able to reduce the impact both of the even zonal 
harmonics of the gravitational field of the central mass and of the non-gravitational perturbations. In general, the 
coefficients m and vj which weigh the various orbital elements in the combinations are a compromise between these 
two distinct needs. 



6.1. The gravito-electric precessions and the solar oblateness 

The main idea is to extend this approach to the Sun-planets system. It must considered that, in this case, the non- 
gravitational perturbations do not play any role 8 , while the N— body classical perturbations of gravitational origin are 
important. Let us start with the gravito-electric precessions writing down the following equations 9 



/ rrSMcrc 
°"obs 



jjoVenus 
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°"obs 
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SZ_2 J 20 
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fz^ ars J. 



20 
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N— body 



'LT MLT: 



'LT 



>LT, 



(15) 



nMars ■ nMars , , 



(16) 



^Mcrc 



vMars j i -Mars..^. 
r .2 ^20+7T GE I/QE, 



•Merc t i 'Merc,^ 
r.2 J 20 + 7T G E ^GE, 



(17) 



where dCl^ not , ^obs" 6 * ano - ^^obs* 16 * are the observational residuals of the rates of the nodes, the longitudes of 
the perihelia and the mean anomalies of Mercury, Venus and Mars. It is intended that all kind of data (optical and 
radio) would be used. The residuals 10 should be built up by purposely switching off the solar quadrupole moment and 
the post-Newtonian gravito-electric accelerations in the force models (or leaving in them some default values to be 
subsequently adjusted according to the present strategy) of the orbital processors. Then, the so obtained observational 
residuals would entirely (or partly, if some default values are left in the force models) adsorb just the investigated 
secular effects and other post-Newtonian short-periodic features, i.e. not averaged over one orbital revolution of the 
planet under consideration. In respect to the latter point, it should be noted that the residuals of the mean anomalies 



Indeed, they are proportional to the area-to-mass ratio S/M of the orbiting probe. The area-to-mass ratio falls off as 1/r 
where r is the radius of the probe-a planet, in this case-assumed spherical. 

9 Note that in the right-hand-sides of eqs. 1161 - 1171 1 also the mismodelled parts of the classical TV— body precessions should 
have been included. Since they are small with respect to the gravito-electric effects of interest, as we will show later, we can 
neglect them in the calculations for setting up our combiantions. This, of course, does not mean that the residuals do not 
account also for them. 

10 Note that here we use expressions like 'observational residuals of a Keplerian orbital element', strictly speaking, in an 
improper sense. The true observables, for which it is possible to build up residuals, are ranges, range rates and angles. The 
Keplerian orbital elements are not directly observable: they can only be calculated from an orbit which has been either numeri- 
cally propagated or fitted to real observational data (the observed orbit) starting from given initial conditions. Here with 'time 
series of the residuals of a Keplerian element' we mean the difference between the time series of that Keplerian element related 
to an observed orbital arc and the time series of that Keplerian element numerically propagated on the same orbital arc (equal 
initial conditions) with the given classical or post-Newtonian parameter of interest switched off in the force models. 
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would account, e.g., also for the indirect effects on the mean motions n through the perturbations in the semimajor 
axes a 



3 / GM 

An = - 2] j—Aa. (18) 

There are no secular perturbations induced on a by the other planets. If the classical short-periodic effects on a would 
be of relatively no importance because they would be included in the force models at the best of their accuracy, this is 
not the case for the post-Newtonian ones. The gravito-electric field induces no secular variations on a, as the classical 
planetary perturbations. The short-term shift on a can be calculated from eq.Q and the Gauss equation for the 
perturbed rate of semimajor axis 

^ = —L=[A R eshxf + A T (l + eco S f)}, (19) 
at nyjl - e 2 

where At is the along-track component of the perturbing acceleration. It amounts to 

Aa GE = f Me [14(cos/ ~ cos/) + 10e(cos 2 / - cos 2 /„)] + 0(e 3 ). (20) 



For Venus their amplitudes are of the order of 100 m; for Mars they amount to 2 km; ea. (|18|l would yield periodic 
variations whose amplitudes would be of the order of 0.2-1 " cy _1 for Venus and Mars, respectively. However, such 
harmonic signatures could be fitted and removed from the data over sufficiently long time spans. Indeed, since we 
are interested in the gravito-electric secular trends on 7r and M it should be possible, in principle, to construct the 
residuals by using orbital arcs longer than the sidereal revolution periods of the planets to be used. Then, all the 
high-frequency perturbations would not affect their time series which should, instead, be characterized by the secular 
parts of those Newtonian and post-Newtonian features present in the real data but (partly) absent in the force models 
of the equations of motion in the orbital processors. However, in regard to the possibility of constructing accurate 
time series of observational residuals many years long the following observations must be kept in mind (Standish 
2002). The planetary motions are perturbed by the presence of many asteroids whose masses are quite poorly known. 
Furthermore, it is not possible to solve for the asteroid masses, other than for the biggest few, because there are too 
many of them for the data to support such an effort. As a result, the ephemerides of the inner planets, especially Mars, 
will deteriorate over time; the ephemerides have uncertainties at the 1-2 km level over the span of the observations 
and growing at the rate of a few km/decade outside that span. On the other hand, it must also be noted that the 
sidereal orbital periods of Mercury and Mars amount to 87.969 and 686.980 days, respectively. 

We can consider eas. <|15[) - (|17H as three systems of algebraic linear equations in the three unknowns 11 J20, Mge and 
i/QE- Their solutions can be written as 



sn 



obs 



ore 



ci^vc.™ 3 + c 2 6n™f s = J 20 (n^ CIC + c x fi^ enus + caO 1 ^ 8 ) 
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c x = -M™ ais /M V 2 CDUS = -0.07475, ( 22 ) 

I c[ = -7r^ ar 77T^ crc = -0.00775, 
and 

' f2^ crc + ciJ7Y, onus + C2 r2^ ars = -32808.8816 " cy" 1 , 
Altars + c[M^ UB = -2.1007 " cy-\ (23) 

^GeT + c 'l'^GE rC = L0176 " C y _1 
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This approach is analogous to that employed in the LAGEOS-LAGEOS II Lense-Thirring experiment in the gravitational 



field of Earth (Ciufolini 1996; Ciufolini et al. 1998) 
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The first equation of ea. (|21|) comes from ea. l|15|) solved for J 2 ©; it allows to obtain J 2 © independently of the post- 
Newtonian Lense-Thirring and classical N— body secular precessions which would represent the major sources of 
systematic errors. The second equation of ea. <(21() comes from ea. <|16[) solved for /iqe; it cancels out the secular 
precessions due to J20. The third equation in eq.J2U comes from ea. ((T7|l solved for ^ge; it cancels out the secular 
precessions due to Jiq- In regard to the impact of the residual mismodelled classical N— body precessions on the 
second and the third combinations of ea. 1)21(1 . they should not induce a systematic error larger than the observational 
one (see below) because the expected values of /xge and ^ge are of the order of unity 12 , contrary to J20 which should 
be of the order of 10 -7 . The adimensional parameters J20, /xge and ^ge are estimated by fitting the time series of 
the left-hand-sides of ea. i|21|) with straight lines, measuring their slopes, in " cy -1 , and, then, by dividing them by the 
the quantities of ea. (|23|l which have the dimensions of " cy -1 . Note that the solar quadrupole mass moment would 
not be affected by the indirect effects on n because only the nodes would be used in its determination. Finally, from 
the so obtained values of /^ge and i^gej which are 1 in GTR and in Newtonian mechanics, it is possible to measure 
7 and independently of the solar oblateness and also of each other as 

f 7 = ^(a^ge + ^ge) - f, 

{ (24) 
1/3 = |mge-|^ge + |. 

According to the results of Table 01 it is possible to yield an estimate of the (formal) uncertainty in J20, /zge and vqe 
as 

' ctj 20 ~ 4.5 x 1(T 9 , 

< a, lGE ~ 6 x 10- 5 , (25) 
k ov GE - 3 x 1(T 5 . 

The evaluation for <r MGE accounts also for the fact that ea. l|18|) and the values of a a of Table VI of (Pitjeva 2001a) 
yield a (formal) uncertainty of 3 x 10 -5 . The uncertainties in 7 and (3 would amount to 

( cr 7 - 6 X 10~ 5 , 

(26) 

[ 073 ~ 1 x 10- 4 . 

However, discretion is advised on evaluating the reliability of these results because they refer to the formal, standard 
statistical errors; realistic errors may be also one order of magnitude larger. 

For the most recent determinations of 7 and (3 we have that, according to the frequency shift of radio photons 
to and from the Cassini spacecraft, a 1 = 2 x 10~ 5 (Bertotti et al. 2003). This result, combined with that for the 
Strong Equivalence Principle violating parameter <t v = 4.5 x 10~ 4 from the most recent analysis of LLR data, yields 
a/3 = 1 x 10~ 4 (Turyshev et al. 2003). 

6.2. The gravito-magnetic precessons 

Let us, now, turn our attention to the Lense-Thirring effect and try to set up some combinations with the nodes of 
the inner planets which cancel out the impact of the Sun's oblateness and of the N— body precessions. 

6.2.1. A J20 — J40 free combination 

Let us write down the expressions of the residuals of the nodes of Mercury, Venus and Mars explicitly in terms of the 
mismodelled secular precessions induced by the quadrupolar and octupolar (Rozelot et al. 2004) mass moments of the 

12 Let us quantitatively discuss this point. By using the results for the observed centennial rates tl and lo released at 
http://ssd.jpl.nasa.gov/elem_planets.html in, say, the left-hand-side of the third combination of ea. 121(1 it is possible 
to obtain for it a nominal TV— body shift of 539.6036 " cy _1 . The uncertainty in the TV— body precessions lies mainly in the Gm 
of the perturbing planets, among which Jupiter plays the major role. Now, the relative uncertainty in Jupiter's Gm is of the 
order of 10 -8 (Jacobson 2003); then, a reasonable estimate of the order of magnitude of the mismodelled part of the JV— body 
shift should be 1 x 10~ 5 " cy" 1 . This fi gure must be divided by 1.0176 " cy 1 yielding a relative error in uge of the order of 
10 -5 . Cfr. with the third equations of eg. 1251 . 
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Sun and of the Lense-Thirring secular precessions, assumed as a totally unmodelled feature. It is accounted for by a 
scaling parameter /l«lt which is zero in Newtonian mechanics and 1 in GTR 

' <^ot rcury = ^ orcury <w 2 © + f^ crcury <5j 40 + sif T crcur V T + A Mercur y 

< Stilus = ^Vcnusj j 2q + ftVenus^ j 4q + fiVeixus^ + A Vcnu Sj (27) 
^ JfjMars = ^Mar S(5j2Q + fiMaWJJ^ + Q^f^ LT + A Mars . 

The coefficients Ctj for the first two even zonal harmonics are 

2 



fi 4 = CI 



1(f) 2 ^(Tsin 2 i-4) 



(28) 



i?0 is the Sun's equatorial radius. The quantities A in eq.J23 refer to the other unmodelled or mismodelled effect 
which affect the temporal evolution of the nodes of the considered planets. In the present case they would be the 
N— body gravitational interactions— of Newtonian origin— with the other planets and asteroids of the Solar System. 
From the nominal values of f^N-body quoted in Table|3|it can be noted that their secular components are very important 
perturbing effects which should be known with very high accuracy in order to prevent a fatal aliasing effect on the 
proposed measurement. 

If we solve eq. (127(1 for the Lense-Thirring parameter /zlt it is possible to obtain 

^ob° s rCUry + /c i < ^obs nus + k 2 Sn^ h f s = X L tMlt + (other secular mismodelled effects A), 



(SM ars q Mercury AMcrcuryAMars 

ki = L^m,,, Lav"! = -48.008308, 



cury^Ve„u S _^ Venus 



(29) 



510.404066, 



Z f!V«.u. n M.,._ f! M.„ f! V,n B! 

X LT = £7LT rCUry + fci^LT nus + k 2 n^ rs = 0.002069 " cy" 1 , 

where the numerical values of the coefficients fci and k 2 and the slope -Xlt of the gravito-magnetic trend come from 
the values of Table |3| The meaning of eq. (|29|l is the following. Let us construct the time series of the residuals of the 
nodes of Mercury, Venus and Mars by using real observational data and the full dynamical models in which GTR is 
purposely set equal to zero, e.g. by using a very large value of c. We expect that, over a multidecadal observational 
time span, the so combined residuals will fully show the GTR signature and partly the mismodelled N— body effect 13 in 
terms of a linear trend. The measured slope, divided by -Xlt, yields /ilt which should be equal to one if GTR was 
correct and if the bias from the residual TV— body effect was sufficiently small. If we use the results of Table0]in ca. (|29H 
in order to get an idea of what could be the observational accuracy of such a measurement we get a 19% error. It is 
important to note that such an estimate does not include the systematic error induced by mismodelled secular effect 
due to the classical N— body interactions. 



6.2.2. A J20 — (N— body) free combination 

Ea. H29|) is designed in order to cancel out the contributions of the first two even zonal harmonics of the multipolar 
expansion of the Sun's gravitational potential. However, it is affected by the N— body classical secular precessions 
which, as it turns out from the figures in Table |3 are quite large; their mismodelled part, which would be accounted 
for by the orbital residuals, could be turn out to be very insidious with respect to our proposed measurement of the 
Lense-Thirring effect. 

In view of these considerations, it would be better to set up another observable which is insensitive just to such 
huge classical secular effects. Moreover, from the results of Table [31 and from the evaluations of (Rozelot et al. 2004) 
according to which the possible magnitude of J40 would span the range 10~ 7 - 10~ 9 , the secular precession induced 
by the octupolar mass moment of the Sun is negligible with respect the Lense-Thirring rates. These facts lead us to 

13 Possible aliasing time-dependent iV— body residual effects with the periodicities of the outer planets, mainly Jupiter, should 
average out over a sufficiently long multidecadal time span. However, it would be possible to fit and remove them from the time 
series. 



12 



L. Iorio: The solar gravito-magnetic field and its possible measurement 



design a three-node combination which cancels out the effects of J20 and of the classical N— body precessions and is 
affected by J^q. It can be obtained by substituting Cl 4 with fl c i ass in ea. i|29|) . It is 

f ^o M b c ; cury + + ~ >itmlt, 

-10.441702, 

(30) 

9.765758, 

^ F LT = Q^ erCUry + fc^Venus + fc^Mars = _ .000351 " Cy" 1 , 

The systematic error affecting ea. (13011 should be totally negligible because it would be due only to the higher degree 
multipole mass moments of the Sun. The observational error, according to the results of Table would amount to 
36%. 

7. Conclusions 

In this chapter we have proposed a relatively cheap approach in testing post-Newtonian gravity in the Solar System 
scenario with particular emphasis on some secular gravito-electromagnetic effects on the orbital motion of the planets. 
It is based on a rcanalysis and a suitable combination of the ranging data to planets. As a by-product, also the 
important solar quadrupole mass moment J2© could be accurately measured. 

First, we have explicitly worked out the post-Newtonian gravito-electric secular rate of the mean anomaly of a 
test particle freely orbiting a spherically symmetric central object. Moreover, we have outlined a possible strategy 
for determining simultaneously and independently of each other, the solar quadrupole mass moment J20 and two 
parameters z^qe and /xqe which account for the post-Newtonian gravito-electric secular shifts of the perihelion and 
the mean anomaly, respectively, of planets. They are in Newtonian mechanics and 1 in the General Theory of 
Relativity. They could be expressed in terms of the standard PPN 7 and parameters which could, then, be determined 
independently of each other as well. The usual approach employed in the ephemerides data reductions tests post- 
Newtonian gravity theories as a whole straightforwardly in terms of the PPN parameters involving the simultaneous 
fit of many more or less correlated astrodynamical parameters among which there are also 7, (3 and J2©- In this 
case, instead, we propose to analyze the time series of three suitably linearly combined residuals of f2, 7T and M. of 
Mercury, Venus and Mars built up in order to single out just certain selected Newtonian and post-Newtonian orbital 
effects. This can be done by setting purposely equal to zero (or to some default values to be subsequently adjusted 
with the proposed strategy) the orbital effects of interest in the force models of the equations of motion. By suitably 
choosing the length of the orbital arcs it would be possible to account for the secular terms only. The coefficients of the 
combinations would make each of them sensitive just to one orbital effect, independently of the other ones. By fitting 
the experimental residual signals with straight lines, measuring their slopes in " cy _1 and suitably normalizing them 
would yield the values of J20, /iqe and i^ge- The obtainable formal accuracy would be of the order of 10~ 4 — 10~ 5 
for 7 and (3 and, more interestingly, 10 -9 for J2©. 

Second, we have extended this strategy to the gravito-magnetic Lense-Thirring effect induced by the Sun's angular 
momentum. We have designed two linear combinations which would allow to measure such tiny effect, by construction, 
independently of the aliasing impact of the Sun's oblateness and of the ./V— body precessions. The obtainable accuracy 
would be of the order of 20-30%. 

However, it must be pointed out that many of the proposed tests involve the analysis of the orbital motion of 
Mercury. The quality and the accuracy of its ephemerides should increase in future thanks to the already launched 
Messenger and planned BepiColombo spacecrafts, so that it is reasonably presumable that the precision of the proposed 
tests will increase as well. 

Acknowledgements. L. Iorio is grateful to L. Guerriero for his support while at Bari. Special thanks to E. M. Standish and E. 
Pitjeva for their help and useful discussions and clarifications, and to W.-T. Ni for the updated reference about ASTROD. 
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